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ABSTRACT 

A ordinary differential equation with stochastic parameters, called SPR_SODE model for the spread of dengue 
fever is considered. Critical Values of bifurcation and the boundary of the above said model are discussed. In this paper, 
the different parameters are considered for further analysis. The bifurcation at the characteristic value of the non-linear 

eigen value equation is supercritical if T l > and subcritical if Tj < 0. The equilibrium solution pair in the positive 
octant of R is also discussed. 
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INTRODUCTION 

Dengue is a disease comes under infectious diseases which is in worldwide. The work of a carrier (i.e) the 
medium for transmitting is performed by the mosquito, "Aedes Ageypti [7]. There are so many models for such infectious 
diseases. We need a separate model for such a special disease like dengue fever for better results. The asymptotically 
stable equilibrium points or equilibrium solutions can be defined as the equilibrium solutions in which solutions that start 
"near" them move toward the equilibrium solution [1]. In this work, the SPR_SODE model [2], [3], [4] 
(SPR_Stochastic Ordinary differential Equation model) is considered to analyze further. All the notions of SPR_SODE 
model [2], [3], [4] are taken for further analysis without any change and the same model is given below. 

SPR_SODE MODEL 
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dt 

|[ /F J=4[^L-Mj^i-^([n)[M 
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d_ 
dt 



dt 



34 



S. Dhevarajan, A. Iyemperumal, S. P. Rajagopalan & D. Kalpana 



(A) 



By converting (A) to fractional quantities and denoting each scaled population by small letters, one can 
get, [2], [3], [4] 



d\ ex \ = fh^fmhWtn r Tp l U_ [ex] m [re] 1 



[ex] h +7] h [if\[ex\ 



dt 



?.+\bir] + J\ 



[TP] 



'<) 



[ TP \ =P h + V„ [ TP \ -{V> ID \ + ^DDD\ [TP]\[TP] h -rj h [in h lTP] 



^\if] m = t m [e X ] m -[BIR] m \if] n 



^[TP] m ^Bim m [TP] m -{[DID] m HDDD] m [TP] m )[TPl 



Now, R Q can be defined as, R Q =^ hm ^ mh , [2],[3],[4] 



where and , can be written in mathematical notation as, 

hm mh ' 



9L 



4 + [DID] m + [DDD] m [TP] m fa [TP] m + fa [TP] 



P [TPl* r * n _1 

— P , [DID] m + [DDD] m [TP~] 



4 + VDID\ + [DDD\ \TP\ fa \TP\ + fa \_TP\ 
■ P„„ + (% + lDID] m + [DDD] m [TP]] ) 



T Jfpu [ i?TPi' ^ +% +iDiD]h hddd] " [7p] *r 



(B) 



(Q 



(D) 



In [2,3,4, S. Dhevarajan et.al, 2013], it was proved that the above said model is asymptotically stable. It is also 
proved that the domain, existence and uniqueness of the solution of the above said model. The disease free equilibrium of 
the model is also proved [3]. The solution of the above said model is asymptotically stable [2]. The disease free 
equilibrium solution is also exists [5]. Asymptotically equilibrium points are near to equilibrium points that are near to 
them move toward the equilibrium solution [1]. 

The Existence of Endemic Equilibrium Points of SPR-SODE Model 

The equilibrium equations for (B) are shown below in (E). In this analysis, Here, the terms 
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\. eX \ \}f~\h 'l- rC -l/i '[^L '[^Im '[^"] t^L are re P resentin g their respective equilibrium values and not their 

actual values at a given time t. 



o= p^rJ^Hl [ T p] h-[ex],-[if] h -[re],]- 



^BIR] h + 



_th_ 
[ TP \ 



0= *»[r\- 



[ex] h +7j h [if] h [ex\ (El) 



(E2) 



L h +[BIR\ 



[ TP \\ 



[ rc \ + n h \if] h [TP] h 

0= P h + d h [TP\~ \[DID\ + [TP] h ).[TP] h - % [if] h [TP] h 



o=^M m -[B^] m [?f 
0=[B/tf] m [77>] m -([D/D] m +[DDD] m [r J P] m )[r J P] m 
Define new parameter, = <f> h I <f) m , to obtain 

[tpT o\tp\ I 

' M m ^-[^[l-[^-[r] ft -N]J- 




Ph 



TP 



■[ TP \{ p ^\ +p J rc \]{ l ~^\ -W J-[4 +[fflRU][«L =o 



(E3) 

(E4) 

(E5) 

(E6) 
(E7) 



[ex\+T] h \if\[ex\=o (E8) 



(E9) 



The bifurcation parameter T can be varied, while keeping all other parameters fixed. In terms of the original 
variables, this corresponds to changing <j> h and <p m , while keeping the ratio between them fixed. Consider = (j) h I (j) m . 

[BIR] m -[DID] m 



One can choose the ratio to sweep out the entire parameter space. Hence, [TP] 

Solving for [(/"] in (E6) in terms of [exj m , one can find [j/*] = 

By rewriting the positive equilibrium for [TP] h in terms of ] from (E9) as 



[DDD] K 



(E10) 



(Ell) 



[ TP l=- 



_ [[BIR] h -[DID], -r/ h [if] h ] + ^[BIR] h - [DID] h -r/ h [if] h ] 2 -4[DDD]^ h 



2[DDD] h 



(E12) 
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Using (E12) in (E3), and solving for [Vc]^ in terms of \jf~^ h '■ 



V rc \=- 



(E13) 



(2L h +[[BIR\ +[DID] h -J7 h [if] h ])+^[[BiR] h -v>m\ -n h W\] 2 +WDD\p h 

It is given the nonlinear nature of (E2), which is not possible to solve for \jf~^^ m terms of [^]^ explicitly. 
Now, by using (E12) rewrite (E2), and define the function [^x]^ = g ^Jf~^ h j as 



, v n . h + r) h +\[[[BIR\ +[DID\ -v h W\]] + §[BIR\ -[DID\ -TJ h [if] h J +4lDDD] t p k 

nWJ= ; \}f\ 

It is clear that g (0) = 0. Coin the positive constant g (1) as [^ max ]^- ^ s ^([^]/j) * s a n * ce smootri an d 
continuous function of \jf~\ h with g '([*/] A )> for \jf \ e t°> ^ and \_ ex \ 6 [°' [ e "*™*L]' 111616 6xistS a 
smooth function \jf~^^ = w ^ domain 0, [^-^wJ^J an d ran g e [0, 1]. As g (0) > 0, the smooth function 

^([^l/J wou ^ exten d to some small [^]^ < 0- an ^ can also 1,6 ex P resse d as functions of 

[^]^by substituting ^(fj = y{\_ ex ]i l ) mt0 (El 2) and (E13). Now introduce D l7 the bounded open subset of 



R defined by 

AH 



[ex] 



Ml 1 



VL Jmj 



-lex] <[ex~] <i 

L Am L Jm 



(E14) 



for some [&x]^> and some[gxf] > 0. Define Z as the open and bounded set Z = { T £ R|-M z < T <M Z }. 

This set is defined to include the characteristic values of L, so there is minimum value that M z can have, but M z may be 

(j) (/), 

arbitrarily large with T = -^ ! — L . 

<f> m N, n +AN h 

Determination of Lower Order Terms 

Now it is to determine lower order terms. (E2) can be written as /| [^]/,'[^]/z) = ^> wriere 

(, [ex\- 6 h +i[[[B/R], +[D/D], -<7,[iy],]] + ^[[B»i], -[D/D], -tfjiflj J +1DDD\ Pl [jf]^ 
and use implicit differentiation to write ^(/"J^ = ^([^"J/,) as a Taylor polynomial of the form 
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Where 3 



8 



df 

' eX \ 



df 

in 



6 h +7 lh +\ [[BIR] h +[DID\] + s j[[BIR] h -[DID],] 2 + 4[DDD] hPh 



The first order approximations to the equilibrium equations can be obtained by substituting (E15) into (E13) and 
(E12), and then all three, along with (El 1) and (E10) into the equilibrium equations (E8) and (E9). Hence, 



/l_10 /l_01 
/2_10 /2_01 



\ex\ 

VL Am J 



+ 



l ex \ 

v \ eX \ 

V.L Am J 



\ 2 



Where, /, 



& + % + \ [l BIR l HDID] h ] + ^[[BIR] h -[DID] h f+4[DDD] hPh 



/i.oi = r 



[B/tf] n [£>£>£>],„ 





[[B/tf], - [D/£>]„ ] + *J[[BIR] h - 


[DID\] 2 +A[DDD\p h 




2[Z>D£>]„ 


£l + ^+|[[™] fc+ [^/D] fc ] 


+ <J[[BIR\ - [DID\ f + 4[DDD\ p h 



0, .p» 



L h + \ [t BIR \ + VDID\ ] ] + ^[BIR\ -[DID\f +4[DDD\p h 



f 2 _ 10 = -([BIR] m +Z m ) 

To apply Corollary 1.12 of Rabinowitz [6], we algebraically manipulate (El 6) to produce 
u = TLu + h(T,u) 



Where U = 



l eX \ 

v \ eX \ , 



and L = 



A 



With 



B 





[[BIR\ - [DID\ ] + J[\_BIR] h - [DID] h f + 4[DDD] h p h 




2[DDD] h ([£/*]„+ 4 )[ 


[O h+ ^\[iBIR\ + VDID\^ 


+ ^[BIR\ -[DID\ f +4[DDD\p h 



6 h .Pmh 



(E16) 



(E17a) 



(E17b) 



(E17c) 



(E17d) 



(E18) 



(E19a) 



(E.19b) 



and h(T, u) is 0(u 2 ). The two characteristic values of L by/^ =1/ \/ AB and X 2 =—\l ^/ AB . As both A and B are 
always positive, due to our assumption that [#^] m > [^ ) ^ ) ] m ' A i s rea l an d corresponds to the dominant eigen value 
of L. The right and left eigen vectors corresponding to A l are respectively, 
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V = (>/A -JbJ andw = (sfB sJa} 



(E20) 



For M z > X l , as £ Y, ( X l , 0) e X . We denote the continuum of solution-pairs emanating from ( X l , 0) by 
ITj , where 1~ [j G X- , and from ( X 2 , 0) by E[ 2 , where ri 2 £ X • Let Z t , Z 2 , C/j and t/ 2 are the sets defined by, 

Zj = {r e Z|w such that (T,u) e II,} (E21a) 

C/j ={weF|r such that (r^jell,} (E21b) 

Z 2 = {r e Z\u such that (r,w) e F[ 2 } (E21c) 

C/ 2 = {w e y |r such that (l\w) e IT 2 } (E21d) [Change the variable, Y and Z] 

The part of Y in the positive quadrant of R 2 can be denoted by Y + and defined 
byY + = {( \ eX \h\ eX \ m ) (= Y | [ ex ] h >0aad [ ex ] m >0 }- The mterna l boundary of Y + can be defined by 



dY + = 



\ ex \ 

v \ eX \ , 

V.L Am J 



Y 



[exj ,>0and "| ([ex] ,=0and "| f[ex] =0and 

r i or r i or r i 

lex \ = o lex >0 ex =0 

V. L Am J VL Jm y VL Jm 



Now, it is to expand the terms of the nonlinear eigen value equation (E18) about the bifurcation point, ( X x , 0). 
The expanded variables are 



u = + su m +s 2 u (1> + (E22a) ; T = \ + eT l +sT 2 +... 

h(r,u) = h{A l +£r 1 +£ 2 r 2 +... ,£u m +s 2 u i2) +... ) 

= £ 2 h 2 [\, U V) ) + .... 



Now, consider, 1 x = — 



W.Lv 



(E22b) 
(E22c) 

(E23) 



Now it is to prove that if v and w are the right and left eigenvectors of L corresponding to the characteristic 
value X x , respectively, the bifurcation at T = X x of the nonlinear eigen value equation (El 8) is supercritical if T { > and 

subcritical if T { < 0. 

Evaluating the substitution of the expansions (E22) into the eigen value equation (El 8) at 0(e 2 ), we obtain 
u <2) = X l Lu <2) + Tj Lu (1) + h 2 , which we can rewrite as 



(I - i!L)u <2) = T l Lv + h 2 , 



(E24) 



Where I is the 2 x 2 identity matrix. As X l is a characteristic value of L, (I - X l L) is a singular matrix. Thus, for 
(E24) to have a solution, T l Lv + h 2 must be in the range of (I -X { L); i.e., it must be orthogonal to the null space of the 
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adjoint of (I-zljL). The null space of the adjoint of (I-zljL) is spanned by the left eigenvector of L (corresponding to the 
eigen value w (E20). The Fredholm condition for the solvability of (A. 24) is w. (TjLv + h 2 ) = 0. Solving for 

Tj provides (E23) obviously. If T { is positive, then for small positive e, u > and r > X x , and the bifurcation is 
supercritical. Similarly, if T \ is negative, then for small positive e, u > and T < A l7 and the bifurcation is subcritical. 

Now, it is to prove that For all u G U b [_ex~j h > and [ex~] m > 0. 

From [5], there are no equilibrium points on 3Y + other than \_6X~j h = = 0, so Ul PI 3Y + = 0. From [5], 

close to the bifurcation point {X { , 0), the direction of Ui is equal to v, the right eigenvector corresponding to the 
characteristic value, X 1 . As v contains only positive terms, Ui is entirely contained in Y + . Hence, for all 
u e Ui, \_exj h > and [ex~\ m > 0. 

Equilibrium in the Boundary of the Positive Orthant 

Now it is to prove that the point u = G Y corresponds to X £ R 7 (on the boundary of the positive 

orthant of R ). For every solution-pair ( 1 , u) G 1 1, , there corresponds one equilibrium-pair (1 , X ) G Z x R , where X 
is in the positive orthant of R 7 . 

First it is to show that u = corresponds to X no(iis ■ As [exj^ = [€Xj = 0' From [1] & [6], One can argue that 
the only possible equilibrium point is X nodis . Now it is to show that for every fe Z, there exists an X in the positive 
orthant of R 7 for the corresponding u G U,. (E25) implies [^xfl > and [^] > 0- The equilibrium equation (Ell) 
implies for every positive > there exists a positive ^(/"J . The equilibrium equation for [TP] m has a positive and 

bounded solution, depending only on parameter values (E10). Now, \jf~^^ = y ( \_ ex ]^ ) implies for every positive 
[/''Jh tnere ex i sts a positive ^yj^ . The equilibrium equations (E13) and (E12) show that, for every positive ^Jf^ 
there exists a positive and [TT 1 ]^ . Hence, the point u = e Y corresponds to X nodis e R 7 (on the boundary of the 

positive orthant of R 7 ). For every solution-pair (T, u) e ]\ x , there corresponds one equilibrium-pair (T ,X )e Z x R 7 , 

where X is in the positive orthant of R (E26) 

CONCLUSIONS 

The different parameters of the SPR_SODE model are considered for further analysis. Here it is proved that the 
bifurcation at the characteristic value of the non-linear eigen value equation is supercritical if T 1 > and subcritical if 

r\ <o. [tp\. It is also proved that the point u = e Y corresponds to X ^ G R 7 , on the boundary of the positive 

orthant of R 7 ) and For every solution-pair ( T , u) G Y[ l7 there corresponds one equilibrium-pair ( T , X ) G Z x R 7 , where 

X is in the positive orthant of R 7 . 
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